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INTRODUCTION 


(Received: 27 August 2018; Accepted: 16 January 2019) 


Nemacheilian loaches are a group of freshwater fishes which are distributed in many 
inland waters of Iran. They are complex species from aspects of systematic with many 
species of this family being morphologically undistinguishable. Due to limitations behind 
the morphology-based identification of such species, genetic approaches such as DNA 
barcoding can be helpful to distinguish species, although their genetic relationships are 
not well determined yet. The present investigation provides data on genetic structure of 
some species of Nemachilidae including Paraschistura bampurensis, Oxynoemacheilus kiabii 
and Turcinemacheilus saadi from inland waters of Iran. These species were sequenced for a 
652 base pair region of the mitochondrial cytochrome oxidase subunit I gene. All species 
were identified by cytochrome oxidase subunit I gene sequence. The sequences of P. 
bampurensis showed that specimens collected from two streams (Shapour and Fahlian, 
south-west of Iran), with 1.89% within-species Kimura two parameter distance, shared 
haplotype in neighbor joining tree. The results provided the first genetic evidence for O. 
kiabii, which has recently been recorded as a new species in Iran, and also cytochrome 
oxidase subunit I gene data showed that O. A/abii is completely different from other 
species of Nemacheilidae. The results suggested that cytochrome oxidase subunit I gene 
sequencing can be used to identify species of Nemachilidae from the Iranian streams. 
The picture of Nemachilidae is yet incomplete, and this study is one of the first steps 
towards enhancing genetic understanding of associations among the species of 
Nemachilidae in Iran. There should also be further attempts in genetic or systematic 
research to fill data gaps among fish species in Iran to help establishing a global COI 
barcode database. 


Key words: Cytochrome oxidase, DNA _ barcoding, Nemachilidae, phylogenetic tree, species 
identification. 


Nemacheilian loaches are benthic species that are distributed in many inland waters of Europe, 
northeast Africa and Asia including streams of Iran (Abdoli, 2000; Coad, 2016), however, their 
characteristics are not well described. Although species of Nemacheilidae are not commercially 
important, the family is of high ecological significance and is a main part of river ecosystems in the 
region. In western Asia, the nemacheilid loaches include numerous species that are morphologically 
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similar. Recent findings regarding Nemacheilidae have led to the identification of new species in Iran 
(Coad & Nalbant, 2005; Golzarianpour ef a/, 2011a, 2011b, 2013; Esmaeili e¢ ai, 2014), which 
means there is still considerable ambiguity in terms of identifying species of Nemachilidae. In other 
words, what is known about them, is mostly based on limited data, and previous studies have mainly 
focused on morphometric aspects of the species from this family (Askari & Shabani, 2013; Esmaeili 
et al., 2013; Kolangi-Miandare e¢ a/., 2013). 

Nalbant and Bianco (1998) reported at least 3 genera and 14 species of nemacheilid loaches 
from Iran. However, recent investigations revealed that the species from this family are more diverse 
than the numbers previously known, with evidence reporting over 40 species including 24 loach 
species endemic to Iran (Abdoli ef a/., 2011; Esmaeili e¢ a/., 2014, Freyhof et al, 2014; Mousavi Sabet 
et al., 2015). Despite recent efforts in identifying these fishes, the list is predicted to be incomplete 
and many taxonomic difficulties remain at the species level mainly due to their small size and low 
commercial value (Tang e¢ a/., 2006; Kottelat, 1990). Nonetheless, delineation and identification of 
fish species is a requirement in studies of ecology, natural history and fisheries management. 
Considering the limitations of anatomical and morphological methods of species identifications such 
as the occurrence of a large morphological plasticity or similarity between organisms of the same 
species, the availability of a small amount of biological material and the existence of convergent 
evolution (Pereira ef a/, 2008; Cawthorn ef a/, 2011), using more accurate and cost effective 
approaches is suggested. 

With the development of a new genetic techniques, the field of species identification have 
relayed on the information from the molecular components of cell. Molecular methods and their 
application for genetic conservation and species differentiation studies have been promoted 
remarkably over the past 20 years (Zhang ef a/., 2004; Comi e¢ a/., 2005; Teletchea, 2009; Dettai e¢ a/, 
2011, Chiu e¢ a/, 2012). DNA-based approaches are among the most used methods for molecular 
species identification (Shendure et a/, 2004; Kochzius et a/, 2010; Zhang & Hanner, 2012). The 
method has been suggested to facilitate the process of species discovery as well as to assign 
unknown individuals to species (Hajibabaei e¢ a/., 2007; Herbert & Gregory, 2005, Mohanty ef a/,, 
2015). 

The term DNA barcoding is generally considered as using a short section of DNA from a 
standardized region of the genome for species identification (Hebert ef a/, 2003a; Hebert ef a/, 
2003b; Hebert & Gregory, 2005). The aim of DNA barcoding is producing a database for origins of 
species and also related sequence information of a 648 base pair (bp) segment of mitochondrial 
cytochrome c oxidase I (CO1), and classifying an organism on the basis of a short sequence of its 
mitochondrial DNA (mtDNA) (Hebert ef a/, 2003a, 2003b; Moritz & Cicero, 2004; Wong ef al, 
2011). The mtDNA has many uses in the field of evolution because of its higher mutation rate and 
lower effective population size compared to nuclear DNA (Brown ef a/, 1979). The COI is 
proposed as a standard barcode for animals (Hebert e¢ a/., 2003a) and could identify a large variety of 
species (Steinke et a/. 2005, 2009; Pegg et a/. 2006). DNA barcoding is also a helpful approach to find 
new species of threatened freshwater fishes (Torres ef a/., 2013). Moreover, mt DNA sequences can 
be utilized for studying population genetics and phylogenetics of fish (Peng ef a/., 2004; Liu & Chen, 
2003). 

Therefore, considering difficulties in the identification of Nemacheilidae species based on 
their morphological characteristics, the aim of this study was to provide barcoding data for some 
species of Nemachilidae in order to rapid and accurate species identification as well as looking 
forward to find internal species marker. Hence, this paper attempts to study the genetic situation of 
some Nemacheilidae species, and to prepare DNA barcoding to compare the relationship among 
species. 
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MATERIAL AND METHODS 

Study area and fish collections 

Nemacheilian loaches were sampled from Shapour, Fahlian and Gamasyab Rivers in the west and 
southwest of Iran (Fig. 1). Shapour and Fahlian Rivers are two important permanent rivers in Fars 
province as these rivers are the main water sources in the region, especially for agricultural activities. 
Gamasyab River is one of the longest permanent rivers in Iran. It originates from the northern 
slopes of Garin Mountain called the Garin headwater. It flows through a few provinces (Hamedan, 
Kermanshah, Lorestan and khuzestan) and ends to the Hoor Al-Azim international wetland. Three 
Nemacheilian loaches included Paraschistura bampurensis from Shapour and Fahlian rivers, and 
Oxynoemacheilus Riabii and Turcinemacheilus saad from Gamasyab River (Fig. 2). The fish were 
collected using nets during May to July, 2012. All specimens were fixed in ethanol (96%) and then 
transferred to the lab. 
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FIGURE 1. Map of the rivers (Gamasyab, Fahlian and Shapour) and sampling locations overlaid on 
the digital elevation model (DEM)-derived hillshade (80 m resolution) representing the topography 
of the earth’s surface in the west and southwest of Iran. 


50 IRANIAN JOURNAL OF ANIMAL BIOSYSTEMATICS Vol.15, No.1 


C 


FIGURE 2. Photos showing Paraschistura bampurensis (A; Length: 55 mm), Oocynoemacheilus kiabii (B; 
Length: 53 mm) and Turinemacheilus saadii (C; Length: 57 mm); three species of Nemacheilian 
loaches assessed in this study. 


DNA extraction and PCR amplification 

Five individuals from three species of Nemachilidae were sequenced, and data were downloaded for 
further 23 individuals of Turcinoemacheilae (Teleostei: Nemacheilidae) from NCBI and EMBL 
GenBank database (Table 1). No sequence of O. Aiabii was recognized in GenBanks. Total DNA 
was extracted from pectoral and pelvic fin using the traditional proteinase-K digestion and standard 
phenol/chloroform protocol storing at -20°C (Hillis ef a/, 1996). In order to amplify fragment of 
mitochondrial COI gene, PCR reactions were conducted using primer cocktails of FishF2-5' 
TCGACTAATCATAAAGATATCGGCAC3' and FishR2-5'ACTTCAGGGTGACCGAAGAATC 
AGAA3' (Ward ef a/. 2005). The 25 ul PCR Reaction mixes included 18.75 ul of ultrapure water, 2.25 
ul of 10 x PCR buffer, 1.25 ul ofMgCl: (50 mM), 0.25 ul of each primer (0.01 mM), 0.125 pl of each 
dNTP (0.05mM), 0.625 U of Taq polymerase, and 0.5—2.0 ul of DNA template. Amplifications were 
performed using a Mastercycler® Eppendorf gradient thermal cycler (Brinkmann Instruments, Inc.). 
The thermal regime consisted of an initial step of 2 min at 95 °C followed by 35 cycles of 0.5 min at 
94 °C, 0.5 min at54 °C, and 1 min at 72 °C, followed in turn by 10 min at 72 °C and then held at 4°C. 
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PCR products were visualized on 1.2% agarose gels containing ethidium bromide (10 mg/ml) and 
the most intense products were selected for sequencing. Products were labeled using BigDye® 
Terminator v.3.1 Cycle Sequencing Kit (Applied Biosystems, Inc.) and then sequenced using an ABI 
3730 apillary sequencer following manufacturer’s instructions. As shown in Table 2, 652 bp of COI 
consensus barcodes for each species were treated as discrete units to estimate the pairwise level of 
genetic divergence using K2P correction model (Nei & Kumar, 2000). 


TABLE 1. List of specimens included in the analyses, their accession numbers, origin, and source of 


sequence. 
Taxon Origin Accession Source 
number 
Paraschistura bampurensis. Iran, Shapout River in Fars Province KP342063 This study 
Paraschistura bampurensis. Iran, Shapout River in Fars Province KP342064 This study 
Paraschistura bampurensis. Iran, Shapout River in Fars Province KP342065 This study 
Paraschistura bampurensis. Iran, Shapout River in Fars Province KP342066 This study 
Paraschistura bampurensis. Iran, Shapout River in Fars Province KP342067 This study 
Oxynoemacheilus Riabii Iran, Gamasiab River in Kermanshah Province KP342068 This study 
Oxynoemacheilus Riabii Iran, Gamasiab River in Kermanshah Province = .KP342069 This study 
Oxynoemacheilus Riabii Iran, Gamasiab River in Kermanshah Province | .KP342070 This study 
Oxynoemacheilus Riabii Iran, Gamasiab River in Kermanshah Province | KP342071 This study 
Oxynoemacheilus Riabii Iran, Gamasiab River in Kermanshah Province = .KP342072 This study 
Turcinemacheilus saadii Iran, Gamasiab River in Kermanshah Province = _KP342073 This study 
Turcinemacheilus saadii Iran, Gamasiab River in Kermanshah Province = .KP342074 This study 
Turcinemacheilus Rosswigi Tran Sirvan KJ 179245 Esmaeili e¢ al, 2014 
Turcinemacheilus Rosswigi Itaq Great Zab KJ179255 Esmaeili et a/,, 2014 
Turcinemacheilus Rosswigi Iran Sitvan KJ179258 Esmaeili e¢ al, 2014 
Turcinemacheilus Rosswigi Iraq Little Zab kJ 179260 Esmaeili et a/, 2014 
Turcinemacheilus Rosswigi Iraq Little Zab KJ 179262 Esmaeili e¢ al, 2014 
Turcinemacheilus Rosswigi Iraq Little Zab KJ179265 Esmaeili e¢ al, 2014 
Turcinemacheilus saadiit Iran Karoun KJ179250 Esmaeili e¢ a4, 2014 
Turcinemacheilus saadit Iran Karkheh KJ 179253 Esmaeili et a/., 2014 
Turcinemacheilus saadiit Iran Karoun KJ179257 Esmaeili e¢ ad 2014 
Turcinemacheilus saadit Tran Karoun KJ179248 Esmaeili e¢ af 2014 
Turcinemacheilus saadii Tran Karoun KkJ179261 Esmaeili e¢ al, 2014 
Turcinemacheilus minimus Turkey Euphrates WJ179251 Esmaeili e¢ al, 2014 
Turcinemacheilus minimus Turkey Euphrates KJ179263 Esmaeili e¢ a4, 2014 
Turcinemacheilus minimus Turkey Euphrates kJ 179249 Esmaeili e¢ al, 2014 
Turcinemacheilus minimus Turkey Euphrates KJ 179256 Esmaeili e¢ al, 2014 
Turcinemacheilus hafexi Tran Karoun KkJ179259 Esmaeili e¢ al, 2014 
Turcinemacheilus hafexi Iran Karoun KJ179252 Esmaeili e¢ a4, 2014 
Turcinemacheilus hafexi Tran Karoun KJ 179254 Esmaeili e¢ al, 2014 
Turcinemacheilus hafexi Iran Karoun KJ179264 Esmaeili et a/., 2014 
Paraschistura bampurensis. Iran Baluchestan KJ179269 Esmaeili e¢ a/, 2014 
Paraschistura bampurensis. Iran Baluchestan KJ 179268 Esmaeili e¢ al, 2014 
Paraschistura malapterura. Iran Namak KJ179267 Esmaeili e¢ a/, 2014 
Paraschistura malapterura. Iran Namak KJ179266 Esmaeili e¢ al, 2014 
Turcinemacheilus bahaii Iran Zayandehroud KJ179246 Esmaeili et a/, 2014 
Turcinemacheilus bahaii Iran Zayandehroud KJj179247 Esmaeili e¢ af 2014 
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TABLE 2. Estimates of Pairwise Genetic Distances among Nemacheilidae species under Kimura 2- 
Parameter Model (Kimura, 1980). 


E1_P. bampurensis 

E3_P. bampurensis 0.000 

F.2_P. bampurensis 0.032 0.032 

E4_P. bampurensis 0.032 0.032 0.000 

E5_P. bampurensis 0.032 0.032 0.000 0.000 

F1_O. kiabii 0.172 0172 0.167 0.167 0.167 

F2_O. kiabii 0.172 0172 0.167 0.167 0.167 0.000 

F3_O. kiabii 0.172 0172 0.167 0.167 0.167 0.000 0.000 

F4_O. kiabii 0.172 0.172 0.167 0.167 0.167 0.000 0.000 0.000 

F5_O. kiabii 0.173 0.173 0.169 0.169 0.169 0.002 0.002 0.002 0.002 

G1_T. saadii 0.158 0158 0.151 0.151 0.151 0.190 0.190 0.190 0.190 0.192 
G2_T. saadii 0.158 0158 0.151 0.151 0.151 0.190 0.190 0.190 0.190 0.192 0.002 
Data Analysis 


The sequences were edited for correction with the SeqScape version 2.6 software (Applied 
Biosystems), and then submitted to the GenBank Barcode database with the accession numbers 
KP342063-74. The sequences from GenBank and BOLD databases and also from our dataset were 
aligned using Mega 5.0 ([amura e¢ a/. 2007). Afterwards, pairwise genetic distances were quantified 
according to Kimura 2-parameter (K2P) distance model (Kimura, 1980) for sequence comparisons. 
Neighbor-joining (NJ) trees of K2P distances were created to provide a graphic representation of 
the patterning of divergence among species (Saitou & Nei, 1987). The robustness of the NJ tree was 
assessed as described by Polgar et a/. (2017). 


RESULTS 

The mitochondrial cytochrome oxidase I (COI) region of samples were successfully amplified using 
PCR. The resulting phenogram of 37 sequences were obtained (Fig. 3). The read lengths were all 652 
bp long. According to the Neighbor-joining (NJ) tree (Fig.3), the species in this study were clustered 
independently. Bootstrap values of species separations were mostly around 100 for O. kzabii and P. 
bampurensis, and around 70 for T. saadiz. NJ tree was consistent in defining the separation among O. 
kiabii and other species, whose clusters were supported by the high bootstrap values (Fig. 3). 

Genetic variations were found within the species with a mean K2P distance of 0.006-1.89. Five 
individuals of O. abit could not be separated from each other. The K2P distance between species 
ranged from the least value at 15.36 (P. bampurensis and T. saadii) to a maximum value at 19.06% (O. 
kiabu and T. saadii). 
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FIGURE 3. Phylogenetic consensus tree of Nemacheilidae species constructed with 652 nucleotide 
of cytochrome oxidase I (COT) gene using Neighbor-joining (NJ) Method. Paracobitis malapterura was 
used as an out group. Bootstrap values greater than 50 are shown. 
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DISSCUSSION 

This study has validated the efficacy of COI gene for identifying species of Nemachilidae family. 
Specimens from three species of Nemachilidae were sequenced for the barcode region of COI. With 
no exception, all samples of the sequenced species were recognized. 100 % of all species was 
amplified with DNA barcoding primer. Three specimens of T. saadii failed to amplify, which could 
be due to DNA degradation. The phylogenetic tree reconstruction suggested that Twrcnemacheilus 
minimus is most closely related to Turcinemacheilus Rosswigi from Tigris drainage, and T. minimus is sister 
group of T. kosswige, After it had initially been identified as T. Rosswigi (Breil & Bohlen, 2001), it 
introduced as a new species, T° minimus, from Turkey (Esmaeili e¢ a/. 2014). Although Tureinemacheilus 
bahaii thas been recorded as a new species in Zayanderoud River from Isfahan, Iran and 
Turcinemacheilus hafexi in Bashar River from Yasouj, Iran (Esmaeili e¢ a/, 2014), these species are 
closely related to each other as sister taxa. T. hafezi is distributed in Bashar River which is connected 
to Karoun and Dez Rivers drainages (Golzarianpour e¢ a/, 2013), flowing to the deltaic area of 
Arvand River in Khuzestan province, Iran. 

The individuals of T. saadiz, which were collected from Sepidbarg and Gamasyab Rivers, 
were separated from each other. According to the phylogenetic tree, Sepidbarg River had more 
diversity of T. saadii in comparison with Gamasyab River. Gamasyab and Sepidberg Rivers flow to 
Karkheh and Karoun Rivers, respectively, however, both rivers shared the same haplotype. T. saadi, 
which was from Bashar River (a branch of Karoun River), is a sister group of T. Aosswigi and T. 
minimus. The color pattern of T. saadii is unique among Turcinoemacheilus in the Middle East region, 
and therefore, it is clearly distinguishable from T. Rosswigi and T. minimus (Esmaeili e¢ al, 2014). T. 
Saadi occurs in sympatry with T. hafezi, but these two species are completely separated from each 
other in phylogenetic tree. 

Divergences among individuals of P. bampurensis were detected. The individuals of P. 
bampurensis were collected from two different watersheds including three specimens from Shapour 
River and two specimens from Fahlian River. These specimens of Shapour River were clustered very 
closely with two specimens of the Fahlian River and far from the other species specimens. Deep 
divergence among individuals of a species that had been assigned to single species may be a result of 
the previously unrecognized species or shared haplotype (Ward e¢ a/, 2005). The average within 
species p-distance of P. bampurensis was 1.89 %o which means that these samples might reflect shared 
haplotypes. The average within species p-distance of the P. bampurensis was remarkably larger than T. 
sdadi, so tesults suggested that two populations of T. Saadi and P. bampurensis have a great diversity 
and COI could clearly separate individuals of these populations from each other. 

The Nemeachilus cladogram (Fig. 3) revealed that O. abit and P. bampurensis are more similar, 
however, O. kiabii is a separate clade. The COI data clearly supported that O. &iabii, which has been 
recorded as a new species in Iran (Golzarianpour e¢ a/, 201 1a), is different from other Nemchilous 
species with acceptable bootstrap value (100 %). Indeed, there were great genetic divergence among 
this species and other Nemachilous species. Contrary to P. bampurensis, other conspecific samples came 
from the same area, thus we might have somewhat underestimated the extent of within species 
diversity. It is suggested that sampling should include individuals from different watersheds for 
freshwater fish as previously mentioned by Ward ef a/. (2005). Confusion in taxonomic assignments 
as a result of inter specific hybridization (Verspoor & Hammart, 1991) does not probably occur. 

Here the usefulness of DNA barcoding has been demonstrated as it has previously been 
stated that this approach is a powerful tool to identify marine and freshwater fish species from 
different geographic regions (Hubert e¢ a/., 2008; McCusker e¢ a/., 2013; Mabragana ef a/., 2012; Victor 
et al, 2009; Keskin & Atar, 2013). As we observed, other studies indicated more than 98 % of the 
analyzed species, especially marine species, could be clearly delimitated through DNA barcoding 
(Zhang & Hanner, 2012; Costa ef a/, 2012). Although we observed a high rate of efficiency, in the 
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case of some freshwater fish species, DNA barcoding has been shown less efficient compared to 
marine species (April et a/, 2011). This is likely due to the fragmentation of rivers and lakes from 
continental freshwater networks. Such fragmentation may consequently lead to a more pronounced 
genetic structure among populations and deeper divergence among haplotypes than in the marine 
realm (Ward ef a/., 1994). Nonetheless, this study is one of the first steps towards enhancing genetic 
understanding of relationships among the species of Nemachilidae in Iran, by using a DNA 
barcoding approach. 

In conclusion, this research has successfully identified species of Nemachilidae using some 
novel and robust molecular techniques. The picture of Nemachilidae is yet incomplete. However, 
the results from this study have led us to conclude that COI barcodes could be helpful to identify 
species from the complex family of freshwater fish. The sequences of O. A/abii are first genetic data 
that were submitted to GenBank (KP342068-72). Further genetic or systematic researches should 
also focus on filling data gaps existing among fish species in water bodies of Iran in order to help 
establishing a global COI barcode database. 
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